Integrated Perturbation Theory and One-loop Power Spectra of Biased Tracers 
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Explicit predictions of the integrated perturbation theory (iPT) for power spectra of biased tracers are gener- 
ally derived in the one-loop approximation. The iPT is a general framework of the nonlinear perturbation theory 
of cosmological density fields in presence of bias, redshift-space distortions, and primordial non-Gaussianity. 
Analytic foiTnulas of auto and cross power spectra of biased tracers including effects of redshift-space distortions 
and primordial non-Gaussianity are comprehensively derived and summarized in the one-loop approximation of 
gravitational nonlinearity. Relations to recent work on nonlinear perturbation theories in literature are clarified 
and discussed. 
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I. INTRODUCTION 

The density fluctuations of the universe contain invaluable 
information on cosmology. For example, the history and in- 
gredients of the universe are encoded in detailed patterns of 
the density fluctuations. The large-scale structure (LSS) of the 
universe is one of the most popular ways to probe the density 
fluctuations in the universe. Spatial distributions of galaxies 
and other astronomical objects which can be observed reflect 
the underlying density fluctuations in the universe. 

In cosmology, it is crucial to investigate the spatial distri- 
butions of dark matter, which dominates the mass of the uni- 
verse. Unfortunately, distributions of dark matter cannot be 
directly observed, because the only interaction we know that 
the dark matter surely has is the gravitational interaction. Con- 
sequently, we need to estimate the density fluctuations of the 
universe by means of indirect probes such as galaxies, which 
have electromagnetic interactions. 

Relations between distributions of observable objects and 
those of dark matter are nontrivial. On very large scales where 
the linear theory can be applied, the relations are reasonably 
represented by the linear bias; the density contrasts of dark 
matter 6m and those of observable objects 6x are proportional 
to each other, 6x - bdm, where b is a constant called the bias 
parameter However, nonlinear eff'ects cannot be neglected 
when we extract cosmological information as much as possi- 
ble from observational data of LSS, and bias relations in non- 
linear regime are not as simple as those in linear regime. 

Observations of LSS play an important role in cosmology. 
Shapes of power spectra of galaxies and clusters contain infor- 
mation on the density parameters of cold dark matter Qqdm, 
baryons £?b and neutrinos Qy in the universe. Precision mea- 
surements of baryon acoustic oscillations (BAO) in galaxy 
power spectra or correlation functions can constrain the nature 
of dark energy [1-3], which is a driving force of the acceler- 
ated expansion of the present universe. The non-Gaussianity 
in the primordial density field induces a scale-dependent bias 
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in biased tracers of LSS on very large scales |]4|-|81. Cosmo- 
logical information contained in detailed features in LSS is so 
rich that there are many ongoing and future surveys of LSS, 
such as BOSS ^, FMOS FastSound Ufll, BigBOSS [H, 
LSST 1 12], Subaru PFS iH, DBS H, Euchd tS, etc. 

Elucidating nonlinear effects on observables in LSS has 
crucial importance in the precision cosmology. While 
strongly nonlinear phenomena are difficult to analytically 
quantify, the perturbation theory is useful in understanding 
quasi-nonlinear regime. The traditional perturbation theory 
describes evolutions of mass density field on large scales 
where the density fluctuations are small. However, spatial dis- 
tributions of astronomical objects such as galaxies do not ex- 
actly follow the mass density field, and they are biased tracers 
of mass density field. Formation processes of astronomical 
objects are governed by strongly nonlinear dynamics, which 
cannot be straightforwardly treated by the traditional pertur- 
bation theory. 

Even though the tracers are produced through strongly non- 
linear processes, it is still sensible to apply the perturbation 
theory to study LSS on large scales. For example, the biasing 
effect in linear theory is simply represented by a bias param- 
eter b as described above. However, biasing effects in higher- 
order perturbation theory are not that simple. A popular model 
of the biasing in the context of nonlinear perturbation theory 
is the Eulerian local bias lll6l420ll . This model employs freely 
fitting parameters in every orders of perturbations, and is j ust 
a phenomenological model because the Eulerian bias should 
not be local in reality. 

The integrated perturbation theory (iPT) ll2lll is a frame- 
work of the perturbation theory to predict observable power 
spectra and any higher-order polyspectra (or the correlation 
functions) of generally biased tracers. In addition, the effects 
of redshift-space distortions and primordial non-Gaussianity 
are naturally incorporated. This theory is general enough so 
that any model of bias can be taken into account. Precise 
mechanisms of bias are still not theoretically understood well, 
and are under active investigations. The framework of iPT 
separates the known physics of gravitational effects on spa- 
tial clustering from the unknown physics of complicated bias. 
The unknown physics of bias is packed into "renormalized 



bias functions" cj in the iPT formalism. Once the renormal- 
ized bias functions are modeled for observable tracers, weakly 
nonlinear effects of gravitational evolutions are taken care of 
by the iPT. The iPT is a generalization of a previous formu- 
lation 1I22I425I1 in which only local models of Lagrangian bias 
can be incorporated. 

In recent developments, the model of bias from the halo ap- 
proach has turned out to be quite useful in understanding the 
cosmological structure formations ll26l - l32ll . The halo bias is 
naturally incorporated in the framework of iPT. Predictions of 
iPT combined with the halo model of bias do not contain any 
fitting parameter once the physical mass of halos are specified. 
This property is quite different from other phenomenological 
approaches to combine the perturbation theory and bias mod- 
els. 

A concept of nonlocal Lagrangian bias has recently at- 
tracted attention 033143511 . Extending the halo approach, a sim- 
ple nonlocal model of Lagrangian bias is recently proposed 
03611 for applications to the iPT. Applying this nonlocal model 
of halo bias to evaluating the scale-dependent bias in the pres- 
ence of primordial non-Gaussianity, not only the results of 
peak-background split are reproduced, but also more general 
formula is obtained. In this paper, the usage of this simple 
model of nonlocal halo bias in the framework of iPT is explic- 
itly explained. 

The bias in the framework of iPT does not have to be a halo 
bias. There are many kinds of tracers for LSS, such as var- 
ious types of galaxies, quasars, Ly-o- absorption lines, 21cm 
absorption and emission lines, etc. Once the bias model for 
each kind of objects is given, it is straightforward to calcu- 
late biased power spectra and polyspectra of those tracers in 
the framework of iPT. As described above, it is needless to 
say that detailed mechanisms of bias for those tracers have 
not been fully understood yet. The iPT separates the difficult 
problems of fully nonlinear biasing from gravitational evolu- 
tions in weakly nonlinear regime. 

While the basic formulation of iPT is developed in 
Ref. piy, explicit calculations of the nonlinear power spec- 
tra are not given in that reference. The purpose of this paper is 
to give explicit expressions of the biased power spectra with 
an arbitrary model of nonlocal bias in the one-loop approxi- 
mations, in which leading-order corrections to the nonlinear 
evolutions are included. The expressions are given both in 
real space and in redshift space. Three-dimensional integrals 
in the formal expressions of one-loop power spectra are ex- 
plicitly reduced to one- and two-dimensional integrals, which 
are easy and convenient for numerical integrations. Contri- 
butions from primordial non-Gaussianity are also taken into 
account in the general expressions. Explicit formulas of the 
renormalized bias functions are provided for a simple model 
of nonlocal halo bias. In this way, the most general formulas 
of power spectra of biased objects in the one-loop approxima- 
tion is provided in this paper 

Since the iPT framework is based on the Lagrangian per- 
turbation theory (LPT) pT-M^, a scheme of resummations of 
higher-order perturbations in terms of the Eulerian perturba- 
tion theory (EPT) [43] is naturally considered il22ll . In this pa- 
per, we clarify the relations of the present formula of iPT and 



some previous methods of resummation technique such as the 
renormalized perturbation theory 1,44. ,45.1 . the Gamma expan- 
sions ll4q - l5lll . the Lagrangian resummation theory ll22l - l25ll . 
and the convolution perturbation theory [52]. Some aspects 
for the future developments of iPT are suggested. 

This paper is organized as follows. In Sec. |II] formal ex- 
pressions of power spectra in the framework of iPT with an 
arbitrary model of bias are derived, and renormalized bias 
functions for a nonlocal model of Lagrangian bias in the halo 
approach are summarized. In Sec.|III] explicit formulas of bi- 
ased power spectra, which are the main results of this paper, 
are derived and presented. Relations to other previous work 
in literature are clarified in Sec. |IV] and conclusions are given 
in Sec. |V] In App. |A] diagrammatic rules of iPT used in this 
paper are briefly summarized. 



II. THE ONE-LOOP POWER SPECTRA IN THE 
INTEGRATED PERTURBATION THEORY 

In this first section, the formalism of iPT is briefly reviewed 
without proofs, and formal expressions of power spectra in the 
one-loop approximation are derived. 



A. Fundamental equations of the integrated perturbation 
theory 

In evaluating the power spectra in iPT, the concept of multi- 
point propagator [46, 47,^3] is useful. The («+ l)-point prop- 
agator ly of any biased objects, which are labeled by X in 
general, is defined by 1I21I1 



S'Sxik) 



66L(ki)---66Ukn) 



= (27T)'-'"6l(k-ki...„)r^^\ku. . . , k„), 

(1) 

where 6xik) is the Fourier transform of the number den- 
sity contrast of biased objects in Eulerian space, d^ik) is the 
Fourier transform of linear density contrast, 6^ is the Dirac's 
delta function in three-dimensions, and we adopt a notation 



ki...„ = ki+' 



l^m 



(2) 



throughout this paper The left-hand side of Eq. ([T) is an en- 
semble average of nth-order functional derivative. The num- 
ber density field is considered as a functional of the initial 
density field. In the basic framework of iPT, the biased ob- 
jects can be any astronomical objects which are observed as 
tracers of the underlying density field in the universe. 

The method how to evaluate multi-point propagators of bi- 
ased objects in the framework of iPT is detailed in Ref. 1I21I1 . 
In the most general form of iPT formalism, both Eulerian 
and Lagrangian pictures of dynamical evolutions can be dealt 
with, and both pictures give equivalent predictions for observ- 
ables. The models of halo bias fall into the category of La- 
grangian bias, i.e., the number density field of halos is related 
to the mass density field in Lagrangian space. In such a case, 
the Lagrangian picture is a natural way to describe evolutions 
of halo number density field. In the models of Lagrangian 



bias, the renormalized bias functions 112111 are the key elements 
in iPT, which are defined by 
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6"6\^{k) 



d^k 



(27r)3 \66dki)---6ddk„) 



(3) 

where 6]^{k) is the Fourier transform of halo number density 
contrast in Lagrangian space. We allow the bias to be nonlocal 
in Lagrangian space. In fact, the halo bias is not purely local 
even in Lagrangian space |36]. Assuming statistical homo- 
geneity in Lagrangian space, the renormalized bias functions 
in Eq. ([3]) is equivalently defined by 1136111 
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FIG. 1 : The diagrammatic representation of the two-point propagator 
with partially resummed vertex up to one-loop contributions. 



d"6\^{k) 



66dki)---66dk„) 



= i2n)'-'"6l(k - ki...„)cf{ku. . . , k„). 

(4) 

The similarity of this equation with Eq. ([U is apparent in this 
form. The information on dynamics of bias in Lagrangian 
space is encoded in the set of renormalized bias functions. 
Assuming statistical isotropy in Lagrangian space, the renor- 
malized bias functions c^ (k\, . . . , k„) depend only on magni- 
tudes ki, . . . ,k„ and relative angles ki ■ kj (i > j) of wavevec- 
tors. 

Applying the vertex resummation of iPT, the multi-point 
propagators of biased objects X are given by a form. 



<«)/ 



r'^\ku...,k„)^n(ki...„)r'^'(k 



,k„). 



(5) 



where 



n{k) = (e-*"^) = exp ^ ^^ {(k ■ W)"\ 

Lh=2 



(-0" 



(6) 



is the vertex resummation factor in terms of the displacement 
field y, and (■ ■ ■ )c indicates the connected part of ensemble 
average. The displacement fields ^(q) are the fundamental 
variables in LPT, where q is the Lagrangian coordinates and 
the Eulerian coordinates are given by j: = q + *P(q). The 
cumulant expansion theorem is used in the second equality 
of Eq. (|6]l. Cumulants of the displacement fields with odd 
number vanish from the parity symmetry, thus the summation 

in the exponent of Eq. (|6]l is actually taken over n = 2, 4, 6, 

The normalized multi -point propagators of the biased objects, 
/J , are naturally predicted in the framework of iPT. 

In the one-loop approximation of iPT, the vertex resumma- 
tion factor is given by 



77^ = exp {-i J 1^ [^ . L(' V)f Pl(p) 



(7) 



and the normalized two-point propagator is given by 

rf{k) = cf{k) + k ■ L^^\k) 

^ ^ ^,Pdp)[^f{k,p)[k-L^\-p)\ 

+ cf{p)[k-L^\-p)\[k-L^'\k)\ 

+ ^kL'^'\k,p,-p) 

+ 4\p)[k-L^^\k,-p)] 

+ [k-L^'\p)][k-L^^\k,-p)]}, (8) 



where L*"^ is the nth-order displacement kernel in LPT. 

Each term in Eq. dHJ respectively corresponds to each di- 
agram of Fig. [1] in the same order Diagrammatic rules in 
iPT |21] with the Lagrangian picture, which are explained in 
App. |A] are applied in the correspondence. The normalized 
two-point propagator of mass density field, tU, is obtained 
by putting c^' = in Eq. (O. 

The perturbative expansion of the displacement field in 
Fourier space, ff'CA:) is given by 






L<-"\ku...,k„)6aki)---6ak„), (9) 



where we adopt a notation. 



(10) 
This notation of Eq. ( fTOb is commonly used throughout this 
paper 

In real space, the kernels of LPT in the standard theory of 



gravity (in the Newtonian limit) are given by ||4] 



r(2)^ 



3 k 



12 






Iky^ 



k\ ■ k2 

k\k2 



L'-'\kuk2, A:3) = ^ [L^'''Hkuk2, ks) + perm.] ; 



(11) 
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Cl23 



k\k2 
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(14) 



where a vector function T represents a transverse part whose 
expression is not used in this paper Complete expressions of 
the displacement kernels of LPT up to 4th order, including 
transverse parts are given in, e.g., Ref. ||4U 14211 . Eqs. (|7) and 
(|8j remain valid even when the non-standard theory of gravity 
is assumed as long as the appropriate form of kernels L„ in 
such a theory is used. 

One of the benefits in the Lagrangian picture is that 
redshift-space distortions are relatively easy to be incorpo- 
rated in the theory. A displacement kernel in redshift space 
U^"'* is simply related to the kernel in real space at the same 
order by a linear mapping ri22|] . 



L<«' ^ L^<"' = L("> +nf[z- L(">) z. 



(15) 



where / = c/lnD/c/lna = D/HD is the linear growth 
rate, D(t) is the linear growth factor, a{t) is the scale fac- 
tor, and H{t) - a/a is the time-dependent Hubble parame- 
ter. The distant-observer approximation is assumed in redshift 
space, and the unit vector z denotes the line-of-sight direc- 
tion. Strictly speaking, the mapping of Eq. ( fTsT i is exact only 
in the Einstein-de Sitter universe. However, this mapping is 
a good approximation in general cosmology. The expressions 
of Eqs. dTji and dS) apply as well in redshift space when the 
displacement kernels in redshift space L^'"' are used instead 
of the real-space counterparts L'"*. 

The three-point propagator at the tree-level approximation 
in iPT is given by 

rf{kuk2) = cf{kuk2) + 4\ki)[k ■ L^'\k2)] 

+ 4\k2) [k ■ L(i)(A:i)] + [k ■ L(i)(A:i)] [k ■ L^'\k2)] 

+ k-L^^\ki,k2), (16) 

where each term respectively corresponds to each diagram of 
Fig. |2]in the same order When the mapping of Eq. (fTST i is 
applied to every displacement kernels in Eq. ( fT6] l, the expres- 
sion of three-point propagator in redshift space is obtained. 
The three-point propagator of mass, /",*„, is given by just sub- 
stituting c^' = in Eq. (fT6T l. 
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FIG. 2: The diagrammatic representation of the three-point propaga- 
tor with partially resummed vertex at the tree-level contribution. 
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FIG. 3: The diagrammatic representation of the power spectrum up 
to one-loop approximation. 



In terms of the multi-point propagators, the power spectrum 
of biased objects, up to the one-loop approximation, is given 
by 



Pxik) ^ n\k)Ur^;^\k)\ Pdk) 

+ i r \r'^^\kuk2)f Pdkl)PL(k2) 



+ r 



f{k) r rf{kuk2)Bdk,kuk2)), 

Jk, 



Jki2=k 



(17) 



where Phik) and B\JJ<.,k\,k2) are the linear power spectrum 
and the linear bispectrum, respectively. The diagrammatic 
representations of Eq. (fTTT i are shown in Fig. |3] Crossed cir- 
cles correspond to the linear power spectrum or the linear bis- 
pectrum, depending on number of lines attached to them. The 
first two terms in Eq. ( fTTj i corresponds to the first two dia- 
grams in Fig.|3] The last two diagrams in Fig.|3]are contribu- 
tions from the primordial non-Gaussianity. The two diagrams 
give the same contribution because of the parity symmetry, 
and the sum of the two diagrams corresponds to the last term 
in Eq. ( fTTT l. 

The matter power spectrum Pmik) is simply given by re- 
placing r^ by r^ in Eq. ( [TtI i. or equivalently, setting c^'^ = 
for every n > 1 . The cross power spectrum between two types 
of objects, X and Y, is similarly obtained as 



PxY(k) = n\k)m\k)r^^\k)pak) 

+ i r rf{kuk2)r^y\k,,k2)Pi.(k{}Pi.(k2) 

^ Jkr = 



-=k 



^ Jk,2=k 

^ '^r?ik) f 
^ Jk,2=k 



r'y'\kuk2)Bdk,ki,k2) 



rf(kuk2)Bdk,kuk2)\. (18) 



The diagrams for the above equations are similar to the ones in 
Fig. [3] where the left and right multi-point propagators corre- 
spond to those of Z and Y, respectively. When X = Y, Eq. (fTsT l 
apparently reduces to Eq. (fTTt . 

The predictions of biased power spectra in the one-loop ap- 
proximation of iPT are given by Eq. (fTTt for the auto power 
spectrum, and by Eq. (fTST i for the cross power spectrum. Once 
a model of the renormalized bias functions c^ is given, it is 
straightforward to numerically evaluate those equations. The 
above results are general and do not depend on bias models. 
Any bias model can be incorporated in the expression of iPT 
through the renormalized bias functions. In the next subsec- 
tion, we explain a simple model of the renormalized bias func- 
tion based on the halo approach. 



in Fourier space. The variance of density fluctuations on the 
mass scale M is defined by 



o-m' 



r cPk 
"J (M 



■WHkR)P^(k). 






j=o J- 



(23) 



The radius R is considered as a function of cr;^i in Eq. ( fT9l ). 
The functions An(M) are defined by 



(24) 



where /jJ; is the scale-independent Lagrangian bias parameter 
of nth-order In this paper, we use first three functions. 



B. Renormalized bias functions in a simple model of halo 
approach 

The renormalized bias functions c^ are arbitrary in the 
general framework of iPT Precise modeling of bias is not a 
trivial problem, depending on what kind of biased tracers are 
considered. In this subsection, we consider the halo bias as 
an example. The expressions of renormalized bias functions 
in a simple model of the halo approach are recently derived 
in Ref. 13611 . We summarize the consequences of this model 
below. It should be emphasized that the general framework of 
iPT does not depend on this specific model of bias. 

Without resorting to approximations such as the peak- 
background split, the halo bias is shown to be nonlocal even in 
Lagrangian space. As a result, the renormalized bias functions 
have nontrivial scale dependencies. For the halos of mass M, 
the renormalized bias functions are given by [i36il 



4V^., 



k„)^^^W{hR)■■■W{knR) 



(5c" 
A„_i(M)crM" d 



Sc" 



d In ctm 



W(kiR) ■ ■ ■ W{k„R) 



o-m" 



(19) 



In this formula, 6^ = 3(3;r/2)^^^/5 - 1 .686 is the critical over- 
density for spherical collapse, W{kR) is a window function 
and the Lagrangian radius R is related to the mass of halos by 



M = -npoR' 



(20) 



where po is the mean density of mass at the present time, or 



R = 



M 



iQi'fmO 



1/3 



/i^'Mpc, (21) 



where M© = L989 x 10 kg is the solai" mass, Q^o is 
the density parameter of mass at the present time, and h - 
//o/(100 km s"' Mpc"') is the normalized Hubble constant. In 
the halo model, the window function is usually chosen to be a 
top-hat one in configuration space, and is given by 



W(x) 



3 sin X - 3xcosx 



(22) 



Ao = l, Ai^l+6,b\, 



A2^2 + 26cb\ + 6^b\. 



(25) 



When the halo mass function n(M) takes a universal form 



n{M)dM = ^/mf(v) — , 
M V 



(26) 



where v - 6c /(Tm, the Lagrangian bias parameters are given 
by 

b];iM) = [^)" 1^, (27) 

where /^'^ = fl'«/MF/«'v". 

Once the model of the mass function, /mf(>') is given, the 
scale-independent bias parameters b]^{M) and the functions 
A„(M) are uniquely given by Eqs. (IZTl i and (l24l i. In Table |I] 
those functions are summarized for popular models of mass 
function, i.e., the Press-Schechter (PS) mass function [26], the 
Sheth-Tormen (ST) mass function [30], Warren et al. (W-h) 
mass function Il54tl . 

In the simplest PS mass function, general expressions of 
the parameters for all orders can be derived [1361 : b^ - 
v"'^Hn+i{v)/6c", An - v"Hn{v), whcrc H„{v) are the Hermite 
polynomials. It is interesting to see that A„ = v^^c" '^J^_i for 
PS mass function. In this case, the renormalized bias function 
c); does not depend on the scale-independent bias parameter 
of the same order, b\, and only depends on lower-order pa- 
rameters b^ , and b^ ^. 

The ST mass function gives a better fit to numerical sim- 
ulations of halos in cold-dark-matter type cosmologies with 
Gaussian initial conditions. The values of parameters in Ta- 
bleUare p = 0.3, q = 0.707, and A(/?) = [\ +n-^l^2-Pr{\l2 - 
/?)]"' is the normalization factor. When we put /? = 0, ^ = 1, 
the ST mass function reduces to the PS mass function. 

The W+ mass function is represented by a parameter cr - 
6clv, which is also a function of M, and parameters A, a, b, c 
are fitted to numerical simulations of halos. The same func- 
tional form is applied to the Marenostrum Institut de Ciencies 
de I'Espai (MICE) simulations in Ref. 1551] . allowing the pa- 
rameters redshift-dependent. Their values are given by A{z) = 
0.58(1 +zr°", a(z) = 1.37(1 +z)-°-^\ b(z) = 0.3(1 -hz)-°°^4^ 
ciz) - 1.036(1 -i-z)"°°^'*. When the redshift-dependent param- 
eters are adopted, the W-i- mass function is sometimes referred 
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Parameters 



A(p) = [1 + ;r-'/-2-"r(l/2 - p)]-' 

p = 0.3 

9 = 0.707 



A-0.13 



W+ : MICE : 

A = 0.7234 A(z) = 0.58(1 +zr 

a = 1.625 fl(z)= 1.37(1 +zr'''5 

b = 0.2538 

c= 1.1982 



ft(z) = 0.3(1 +z) 
c(z)= 1.036(1 +z) 



-0.084 



-0.024 



TABLE I: Functions ^^^(M), A„(M) derived from several models of mass function. 



to as the "MICE mass function". In the latter case, the multi- 
plicity function /mf(v) explicitly depends on the redshift, and 
the mass function is no longer universal. 

The Eq. ( fT9b is shown to be equivalent to the following ex- 
pression. 



ct\ku. 



.,k„)^b\{M)W(hR)- 
A„_i(M) d 



dlncTM 



W{k„R) 
[W{hR)---W{knR)]. (28) 



The nonlocal nature of the halo bias in Lagrangian space is 
encoded in the second term in the RHS of the above expres- 
sion, since the simple dependence on the window function of 
the first term appears even in the local bias models through 
the smoothed mass density field. In the large-scale limit, 
ki,k2,...,k„ -^ 0, the second term in the RHS of Eq. ( l28b dis- 
appears and the renormalized bias functions reduce to scale- 
dependent bias parameters, c^^ ^ b];(M). However, the loop 
corrections in the iPT involve integrations over the wavevec- 
tors of the renormalized bias functions, and there is no reason 
to neglect the second term which represents nonlocal nature 
of Lagrangian bias of halos. 

In the expressions of renormalized bias functions, Eqs. (|T9T l 
and ( |28] l, all the halos are assumed to have the same mass, M. 
These expressions apply when the mass range of halos in a 
given sample is sufficiently narrow. When the mass range is 
finitely extended, the expressions should be replaced by [|36ll 



c';'(ku...,k„) 



JdM4>(M)n(M)c%\ki,...,k„) 
jdM(p(M)n{M) 



(29) 



where n{M) is the halo mass function of Eq. (|26l l, 0(M) is a 
selection function of mass. For a simple example, when the 
mass of halos are selected by a finite range [Mi , M2], we have 






c;.; „.(ki,.. 



> i^n) 



CdMn(M)c%\ku...,k„) 



rM- 



UdMn(M) 



(30) 



The number density of halos are usually decreasing function 
of mass, and the integral of Eq. ( |30] | is mostly weighted at the 
low-mass end. Mi . 



III. EXPLICIT FORMULAS 

The auto power spectrum Px{k) of Eq. ( flTl l is a special case 
of the cross power spectrum Pxrik) of Eq. ( fTsT i as the former 
is given by setting X - Y in the latter It is general enough to 
give the formulas for the cross power spectrum below. In the 
following, we decompose Eq. ( fTSb into the following form: 



Pxrik) = n\k) 

X 



Rxrik) + ^QxY(k) + -SxY(k) + -Syxik) 



(31) 



where IJ(k) is given by Eq. (|7]i and 



(32) 
QxY(k)^ [ rf{kuk2)rf{kuk2)Pdkx)PL{k2\ (33) 



lkr=k 



SxY(k) = r' 



Jkn 



=k 



rf\ki,k2)Bdk,kuk2). 



(34) 



Three-dimensional integrals appeared in the above compo- 
nents of Eqs. (l32Ti-(l34li can be reduced to lower-dimensional 
integrals both in real space and in redshift space. Such dimen- 
sional reductions of the integrals are useful for practical calcu- 
lations. The purpose of this section is to give explicit formulas 
for the above components 77, Rxy, Qxy, Sxy in terms of two- 
dimensional integrals at most. The results of this section are 
applicable to any bias models, and do not depend on specific 
forms of renormaUzed bias functions, e.g., those explained in 
SecHm 



A. The power spectra In real space 

In real space, the power spectrum is independent on the di- 
rection of wavevector k, and thus the components above n{k), 
Rxvik), QxYik), Sxrik) are also independent on the direction. 
In this case, the dimensional reductions of the integrals in 
Eqs. ([32]l-([34ll are not difficult, because of the rotational sym- 
metry. The vertex resummation factor n(k) of Eq. (|7]i is given 
by 



nik) = exp 



Un^ 



f 



dpPdp) 



(35) 



On small scales, this factor exponentially suppresses the 
power too much, and such a behavior is not physical. This 
property is a good indicator of which scales the perturbation 
theory should not be applied. However, the resummation of 
the vertex factor is not compulsory in the iPT. When the ver- 
tex factor is not resummed, one can expand the factor as 



n{k) = 1 - 



Un^ 



/ 



dpPdp), 



(36) 



instead of Eq. (1351 ) in the case of one-loop perturbation the- 
ory. In a quasi-linear regime, the resummed vertex factor of 
Eg. |35]) gives better fit to A^-body simulations in real space 

The expression of two-point propagator in Eq. (|8) is 
straightforwardly obtained, substituting the Lagrangian ker- 
nels of Eqs. (fTTI)-(fT4li. Taking the z-axis of p as the di- 
rection of k, integrations by the azimuthal angle are trivial. 
Transforming the rest of integration variables as r = p/k and 
X - p ■ k,we have two equivalent expressions. 



f<'*(/t) = 1 + c<l*(/t) + 



= 1 + 4\k) + 



^3 /-KJ M 

—^ I dr I dx'Rx(k,r,x)Pdk^ 

k^ r ~ 

-^ I dr'Rx(k,r)Pi^(kr), 



r) 
(37) 

(38) 



where 



'Rxik,r,x)^ 



2\2 



5 r\l-x^) 
21 1+r^- 2rx 
3 (!-«)(! -x^) 



7 1 



• 2rx 



[rx + r^c''^\kr)] 



rxcf(k,kr,x), (39) 



and 



'Rxik,r) = 



6-H5r2-H50r'*-21r'' 
252r2 
(1 - r^)\2 + Ir^) 
"^ 1687^ 

3 -H 8r2 - 3r'* 3(1 - r^f 



28 



56r 



1 -r 



1 +r 



■x: 



1 +r 



„(2) 



4^'(fer) 



dxxc'"'ik,kr,x). (40) 



In the above expressions, rotationally invariant arguments for 
c^ are used, i.e.. 



cy (ki,k2) = cy(kuk2;x), 



(41) 



where x = ^i ■ ^2 is the direction cosine between ki and ^2- 
The second expression of Eq. ( |38] | is obtained by analytically 
integrating the variable x in the first expression of Eq. (|37] |. 
Both expressions are suitable for numerical evaluations. With 
the expression of Eq. (l37T i or (l38T l. we have 



RxYik) = r<^\k)r^^\k)Pdk). 



(42) 



Evaluating the convolution integrals in Eqs. ( l33T l and (l34l i 
with the three-point propagator of Eq. (fT&t is also straight- 
forward in real space. Substituting the Lagrangian kernels of 
Eqs. ( fTTT i and ( fT2l i into Eq. ( fT6] l, and transforming the integra- 
tion variables as r = fci/fc, x = ^ ■ ^1, we have 



QxYik) 



^r*x: 



dx p-rfiK r, x)rf{k, r, x) 

X Pdkr)Pdky) (43) 



and 







<i)/ 



Sxrik) = ^2^'x'(k) 

X 



XOO pi 

dr I dxr^rf\k,r,x)Bi^(k,kr,ky) (44) 



where 



and 



Vl + r2 - 2rx, 



(45) 



ff(A:,r,x) 



4 1 
I — rx 



X X 

+ 
r 



-[1 

r L 



.(1) 



(ky) 



[1 + cf(kr)] + cfikr, ky; x), (46) 



The factor ly (k, r, x) is simply given by substituting X ^ Y 
in Eq. ( |46] |. The function t^ (k,r,x) is just the normalized 
three-point propagator P^ {k\,k - ^1) as a function of trans- 
formed variables. 

All the necessary components to calculate the power spec- 
trum of Eq. (ISTT i in real space. 



PxYik) = n\k) 

X 



RxYik) + -QxY(k) + -SxYik) + -SYx(k) 



(47) 



are given above, i.e., Eqs. ^^ [or ^B], 62), dlSll and (Hil l. 
Numerical integrations of Eqs. (OD [or (gT)], (g3) and (g4) are 
not difficult, once the model of renormalized bias functions 



„(«) 



Cj^ and primordial spectra Piik), B]JJ<.\ , k2, k^) are given. The 
last factors SxYik) and Syxik) are absent in the case of Gaus- 
sian initial conditions. 
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T{k,p) 




/(2;^3mp)/'.(p) 




Diagram 


V^\k,p,-p) 


^P-<« 


n^ 


Lf\-p)Lf(k,p) 


3 


kikj - k^Sjj 3 kikj 
k^ ^'(^>-7^^ 


2^ 


lyj-- 


V'^\k,p)4\p) 


7P^3(^) 


[yX- 


L^'\-p)cf(k,p) 




- I'^^(^) 




^ 



TABLE II: Integral formulas for one-loop corrections, which are re- 
lated to the two-point propagator. We denote Ri(k) = R^(k) and 
R2(k) = i^f (fc), as these functions are independent on the bias. 



B. Kernel Integrals 

Evaluations of power spectra in redshift space are more 
tedious than those in real space. The reason is that the 
power spectra depend on the lines-of-sight direction in red- 
shift space. One cannot arbitrary choose the direction of z-axis 
in the three-dimensional integrations of Eqs. dHJ and (l32t- 
( |34] i. because the rotational symmetry is broken. Even in such 
cases, an axial symmetry around the lines of sight remains, 
and the three-dimensional integrations can be reduced to two- 
or one-dimensional integrations as shown below. All the nec- 
essary techniques for such reductions are the same with those 
presented in Refs. ll22ll23ll . making use of rotational covari- 
ance. We summarize useful formulas for the reduction in this 
subsection. We assume the standard theory of gravity in the 
formula below, although the same technique may be applica- 
ble to other theories such as the modified gravity, etc. 

The first set of formulas is related to the two-point propa- 
gator r^ of Eq. (O. The results are summarized in Table llll 
The integrals of a form. 



/ 



(2^ 



nk,p)Pdp), 



(48) 



where T^ik, p) consists of LPT kernels L'"' and renormalized 
bias functions c^ , are reduced to one-dimensional integrals, 
R^{k). The explicit formulas are given in Table HH In this Ta- 
ble, we denote Ri{k) = R^(k) and Riik) = ^f (/fe), as these 
functions are independent on the bias. The functions R^{k) 
are defined by three equivalent sets of equations. 



Rn(k) 



(k,p)Pdp) 



^J drnl{r)P^{kr), 



dx'R^(r,x)Pi^(kr) 



where integrands "R^ik, p), "R^ir, x), and ^^(r) are given in 
Table Unl The last expression of Eq. ( l49b is the formula which 
is practically useful for numerical evaluations. The other ex- 
pressions are shown to indicate origins of the integrals. 

When the second-order bias function c\^ (ki,k2) only de- 
pends on magnitudes of wavevectors ki and k2, and not on the 
relative angle pi2 = ^i ■ ^2, the fourth function genetically 
vanishes: Rf(k) — 0. When the first-order bias function c^ 
is scale-independent, it is explicitly shown from the last ex- 
pressions that R^{k) - [R\{k) + R2{k)]c'^\ Specifically, the 
functions R^{k) and R^{k) are redundant in the Lagrangian lo- 
cal bias models, in which renormalized bias functions cj; are 
scale-independent. This is the reason we only need R\{k) and 
R2(k) in Ref. ||23I1 . In general situations with Lagrangian non- 
local bias models, aU four functions are needed. In a simple 
model of halo bias in this paper, the second-order bias func- 
tion cy does not depend on the angle ^n and Rn{k) = in 
this case. 

The second set of formulas is related to the convolution in- 

(2) 

tegrals of the three-point propagators F^ in calculating the 
one-loop power spectrum. The integrals of the form. 



r 

Jk,- 



T(kuk2)PL(kx)Pak2\ 



(50) 



where T consists of LPT kernels L„ and renormalized bias 
functions c^ and Cy, are reduced to two-dimensional inte- 
grals, QfJ{k). The explicit formulas are given Table HV] For 
the third and fifth formulas in this Table, the indices of the 
LPT kernels are symmetrized, since only symmetric com- 
binations are used in this paper. In this Table, we denote 
Qn{k) - QfJ{k) for n - 1, 2, 3, 4, as these functions are inde- 
pendent on the bias, and g^(/t) = (^/{k) for n = 5, 6, 7, 8, 9, 
as these functions are only dependent on the bias of objects X. 
The functions Q^^ik) are defined by two equivalent sets of 
equations. 



ef(fe)= r Q'/(kuk2)Pdki)PL(k2) 

Jk,2=k 

dr I dxQ^^(r,x)Pi^(kr) 

xPi^[ky/l +r^-2rx), (51) 



Jki2=k 
k^ r- 

4^2 



where integrands Q^^(k[,k2), Q^^{r, x) are given in Table IV] 
The last expression of Eq. ( fsTl i is the formula which is prac- 
tically useful for numerical evaluations. The first expressions 
are shown to indicate sources of the integrands. 

The third set of formulas is related to the initial bispectrum, 
which is an indicator of primordial non-Gaussianity. The in- 
tegrals of the form. 



/ 

Jk,-. 



r(ki,k2)Bdk,kuk2), 



(52) 



where f consists of LPT kernels L„ and renormalized 

bias functions c^ , are reduced to two-dimensional integrals, 

(49) S^(k). The explicit formulas are given in Table fVll In this 

Table, we denote S \(k) = Sf(k) and S2ik) - ^^(A;), as these 



'K(k,p) 



\k-p? 



1- 



{k-p)[k-{k-p)\ 



p-lk-pp 



k_p 

kp 



1- 



k(k-p) 

\k-p? 



k P 

kp 



k ■ p 

~kf 



^'(P) 



^cfik,p) 



rHl - ^f 
1 + r2 - 2rx 

rx(l - rx)(l - x^) 



1 + r2 - 2rx 



rHl - rx)(l - x2) ,„ 



I + r^ — 2rx 



c'i'ikr) 



rxc\/{k,kr,x) 




cf{kr) 



r I dxxc\\k,kr\x) 



TABLE III: Integrands for functions R^ik) of Eq. ( |49t . 



functions are independent on the bias. The functions 5^(fe) 
are defined by two equivalent sets of equations. 



Jk,2=k 
,3 ^oo ^1 



: Bl {k, kr, k Vl + r^ - 2rx) , 



(53) 



where integrands S^(ki,k2), S^{r, x) are given in Table I VII" 



C. The power spectra in redshlft space 

As all the necessary integral formulas are derived in the pre- 
vious subsection, we are ready for write down the explicit for- 
mula of the power spectrum in redshift space. The decompo- 
sition of Eq. ( |3TI ) is applicable in redshift space, and it is suf- 
ficient to give the explicit expressions for the functions n{k), 
Rxrik), Qxrik), S xrik) in redshift space. These functions 
depends on not only the magnitude k but also the direction 
relative to the lines of sight. 

We employ the distant-observer approximation for the 
redshift-space distortions, and the lines of sight are fixed in the 
direction of third axis, z. Lagrangian kernels are replaced ac- 
cording to Eq. ( fTSl l in the formulas of propagators Eqs. (JS) and 
( fTSI l. In those formulas, the Lagrangian kernels appear only in 
the form of k ■ L*"\ With the linear mapping of Eq. ( flSl l. we 
have 



where 



k ■ L^"> -^ k ■ L'^"> = (k + nffikz) ■ L^"\ (54) 



fi = k ■ z, (55) 



and k = k/k. Thus, in the distant-observer approximation of 
this paper, the direction dependence comes into the formulas 
only through the direction cosine of Eq. (ISST l. We denote the 
functions of Eqs. (I32b-(l34l) as RxyiK^i), QxriK^i), SxY(k,^i) 
in the following. 

Substituting Eq. ( |54b into Eqs. Q, ® and ( fT6] l, one can 
see that evaluations of Eqs. (I32]i-(l34li are straightforward by 



means of the integral formulas in the previous subsection. The 
results are explicitly presented in the following. 

The vertex resummation function of Eq. ^ can be evalu- 
ated by applying the same technique of the previous section. 
The relevant integral is 



/ 



^(.-.•■■V)1V,(., 

^ (ki + fiukzi)(kj + ffikzj) I — -^ 



-^^Pl(p), (56) 
y P^ 



and the last integral is proportional to the Kronecker's delta. 
The proportional factor is evaluated by taking contraction of 
the indices. Consequently, we have 

n(k,p) = exp |- [l + /(/ + 2)iJ^] -^ jdpPdp)] ■ (57) 

The two-point propagator of Eq. dH) with the substitution 
of Eq. ( |54| | is evaluated by means of Table [III where R„{k) 
functions are defined by Eq. ( |49] l and Table |III] The result is 
given by 



l + -Ri + -R2 + -Ri-Rf 



//^^ 



7^i/V + 



3 6 

-Ri + -R2 

7 7 



fV 



(58) 



„(i) 



The quantities c^ , R„, R^ on LHS are functions of A:, although 
the arguments are omitted. The component Rxy of Eq. (l32b 
is straightforwardly obtained by the above result of the two- 
point propagator: 



RxY(k,p) = r'^^\k,i^)r^^\k,i^)Pdk). 



(59) 



The tree-level contribution of the above equation is given by 
(bx+ffi^)(bY+ffJ.^)Ph(k)v/hevebx - l-nc^'', and Kaiser's lin- 
ear formula of redshift-space distortions for the power spec- 
trum 1 56] is exactly reproduced. In calculating the mass power 
spectrum, X - Y - m, we only need terms with Ri(k) and 
R2(k) and other terms R^(k) and R^(k) vanish since c*^'* = 
for unbiased mass density field. 
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T{k,,k2) 



r r(k,,k2)Pdki)PL(k2) 



Diagram 



LY\kuk2)Lf{ki,k2) 



^^e.(^) 



I [w\/^ y^^v/v\j I 



L\'\ki)L\''>(k2)Lf'>(kt,k2) 



3 (kikj-l^5n)kk 'ihkjh 



^CJ^"'^ 



3 kjkjkuk, - Ik^S^ijhk,) + k'^SajSu) 



Qi(k) 



L'^(k^)L\'\ki)L['\k2)L]l\k2) 



1 kjkjk^ki — k^6{ijktki-f kikjk^ki 



<l> 



fcs 



-Q3(k)+—T^Q4(k) 



Lf\ki)Lf(ki,k2)c'-^\k2) 



3 kjkj ^^^ 
7 /(^ 



Lj,"(A:,)Ly'(A:i)4';(A:2)4"(fc2) 



1 /C y /C y /C^ /C Oij j K-j^ \ Y '^i '^ j '^k 



k6 



Qi(k) + ^Q^(k) 



°C:> 



V^\k,,k2)cf{k„k2) 



'~QUk^ 



L\'\k,)Lf{k2)cf{kuk2) 



1 KjK j K. j j V- "^i '^ i V 

2 y ' e;m->ye;( 



o 



Lf\k,)L\'\k2)4\ki)c^^\k2) 



1 kikj-k^6i^ 

2 ^4 



efo-'w + ^effw 



o 



L<"(A:,)L<"(Ai)4'>fe)4"(^2) 



\_hk^-je^ 

' 2 k^ 



Qn{k)^^-^Q\l{k) 



nc::^ 



L''\ki)c^'\k2)cf(k,,k2) 



leflw 



<> 



cf{kuk2)cf{kuk2) 



0l(.k) 



o 



TABLE IV: Integral formulas for one-loop corrections, which are related to convolving three-point propagators. In the third and fifth formulas, 
the spatial indices are completely symmetrized. We denote Q„(k) = Q^ik) for n = 1, 2, 3, 4, as these functions are independent on the bias, 
and Q'^ik) = Q^/(k) for n = 5,6,7, 8, 9, as these functions are only dependent on the bias of objects X. 



The component QxY{k,n) of Eq. dSST l is similarly evaluated, 
while the number of terms are larger The result is given by 

QxYiK fD^Yj J"'"/'" [«™ (^) + «™ (^)] (60) 



where 



«TO = ^Gi + 7G2 + -24 + -e^ + 22^ + -Ql + Q; 



X , r\X 



+ efr + ef3''+2e- + ief/, 



(61) 



(62) 



+2eff+2ef/+2efi'. 



^12 = -y^Gi + 463 + G^ - ^Qt - ifiXl + ^Q\h (63) 



.XF_57^ ,15^ 1^ , ,^ , 12^x ^z 



^2 = 98 Gi + y G2 - 4G3 + 324 + y Gf - Ge + 6G7 

(64) 
(65) 



^\q1^qI^\q'^,^q\1-\q\1^q\1 



all = -7G1 + ^G3 + Ge 



16^ 

3 



flf/^TrGi, 



4/ = ^Gi + ^G2 - ^G3 + 224 - Ge + 2Gf , 



^34 



^-^1' 



3 1 

-3G1 + 4G3, 

1 



44-^Gl-4G3 + -G4, 



(66) 
(67) 
(68) 
(69) 



and other a^lSK) which are not listed above all vanish. The 
quantities 2„, Q^, Q^^ are functions of k, although the ar- 
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kiko 



cf(ki,k2) 
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e 



1- 



kki 



kk. 



C''>(k2)c';'(k2) 



kk, 



i M>n..\JM 



4"(^2)d"(fc2) 



kk 



' „(i)^;, \,.(2), 



4'(k2Wy"(kl,k2) 



ki' 

cf(ki,k2)cf(ki,k2) 



r' c\^ (kr,ky;fj.) 
X cf\kr, 

TABLE V: Integrands for functions QY(k) of Eq. dSTT l. 
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Diagram 
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TABLE VI: Integral formulas for one-loop corrections, which are 
related to convolving three-point propagators with the linear bispec- 
trum. We denote S \(k) = S'^ik) and S2(k) = S^C^-), as these func- 
tions are independent on the bias. 



S^(ki,k2) 

[k = k, + k2] 




k,%' 



kk 



-4"(fc2) 



ki' ' 

cf{kuk2) 



S^'Sr^x) 
y = (\ + r^- 2rx)"^ 
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~J' 
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y' 



rxc^ \ky) 
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TABLE VII: Integrands for functions Sl(k) of Eq. | |53> . 



guments are omitted. The Qn functions of n - 1 , . . . , 4, 
10, . . . , 13, 15 are symmetric with respect to X <-^ Y, while 
those of n = 5, . . . , 9, 14 are not. In calculating cross power 
spectra, X i^ Y, the symmetrization with respect to XY in 
Eq. (|60] l is necessary. In calculating auto power spectra, 
X - Y, two terms in the square bracket in Eq. (l6Qt are the 
same, and can be replaced by 2a^^^{k). In calculating the 
mass power spectrum, X = F = m, we only need terms with 
Qiik), ..., Q4(k) and other terms Q^(k), ..., Q^^(k) all vanish 
since c^ = for unbiased mass density field. 

The component 5 xrC^,/") of Eq. ( [34] i is similarly evaluated. 
The result is given by 



SxY{k,ti)^r^^\k,ii) 



-Si +S2+2S^ + S 



+ {-Si+2S2+2Si\fij' 



f2,.2 



--5i/V+ 2^1+521/^ 



2, .4 



(70) 



The quantities S „ and S^ are functions of k, although the ar- 
guments are omitted. The normalized two-point propagator 
t]^ (k,ju) in Eq. ( fTOb can be replaced by the tree-level term, 
1 + 4 + ffi^, because the rest of the terms are orders of two- 
loop approximation. 

All the necessary components to calculate the power spec- 
trum of Eq. (n\ in redshift space. 



PxY(k,fi)^n\k,ti) 

RxYik,n) + -QxY(k,ti) + -SxY(k,fi) + -SYx(k,^i) 



(71) 



are given above, i.e., Eqs. (|57T i, d59] l, ( |60] l and (fTOl i. Numerical 
integrations of Eqs. (|49] l, (ISTT i and ( |53] l are not difficult, once 
the model of renormalized bias functions c^ and primordial 
spectra P^ik), Bi^(ki,k2, kj,) are given. The last factors Sxyik) 
and S Yxik) are absent in the case of Gaussian initial condi- 
tions. 
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D. Evaluating correlation functions 

We have derived full expressions of power spectra of biased 
tracers in one-loop approximations. The correlation func- 
tions are obtained by Fourier transforms of the resulting po wer 
spectrum. In real space, the relation between the correlation 
function ^xrif) and the power spectrum Pxrik) is standard: 



^XY(r) : 



f 

Jo 



k^dk 

1^ 



■Mkr)PxY(k), 



(72) 



where ji(z) is the spherical Bessel function. For a numerical 
evaluation, it is convenient to first tabulate the values of power 
spectrum Pxrik) of Eq. WT\ for a one-dimensional integration 
ofEq. (III). 

In redshift space, multipole expansions of the correlation 
function are useful ll57l - l59ll . For reader's convenience, we 
summarize here the set of equations to numerically evaluate 
the correlation functions in redshift space from the iPT for- 
mulas of power spectra derived above. The multipole expan- 
sion of the power spectrum in redshift space, PxY{k,iJ), with 
respect to the direction cosine relative to lines of sight has a 
form. 



PxY(k,tl)^YjP'xYik)Pl(l^\ 



(73) 



/=() 



where P/(//) is the Legendre polynomial. Inverting the above 
equation by the orthogonal relation of Legendre polynomials, 
the coefficient p'xyik) is given by 



PxY^k) 



21+ 1 



r dfiP,(pi)PxY(k,iJ). 



(74) 



Because of the distant-observer approximation, the index I 
only takes even integers. 

The dependence on the direction fi of our power spectrum, 
PxYik,//) of Eq. ( TTTT i, appears in forms of yU^"e""''' where 
n = 0, 1, 2, . . . are non-negative integers. It is possible to ana- 
lytically reduce the integral of Eq. (l74t by using an identity 



x: 



diui/"e 



2«.-'V ^a-"-^''-y\n + -,a 



(75) 



where y{z, p) is the lower incomplete gamma function defined 
by 



y{z 



,P)^ f'e-'i 
Jo 



t ^dt. 



(76) 



Although the number of terms is large, it is straightforward to 
obtain the analytic expression of p'^yik) of Eq. ( l74l i in terms of 
Q„(k), R„(k), S„(k), c^ (k), and the lower incomplete gamma 
function. Computer algebra like Mathematica should be use- 
ful for that purpose. Alternatively, it is feasible to numer- 
ically integrate the one-dimensional integral of Eq. (l74l for 
each k, once the functions Q„(k), R„ik), S„ik), c^^ (k) are pre- 
computed and tabulated. The latter method is much simpler 
than the former. 



The multipole expansion of the correlation function in red- 
shift space, ^xYir, fj), with respect to the direction cosine rela- 
tive to lines of sight is given by 



aY(r,^)^Yj^'xY(r)Pi(Ml 

1=0 

, 21+1 r^ 

^xy(>-) = -^- J dfiP,(fi)^xY(r,l^)- 



(77) 
(78) 



Since the power spectrum PxY(k,l-i) and the correlation func- 
tion ^xYik, fj) are related by a three-dimensional Fourier trans- 
form, corresponding multipoles are related by [57] 



^XY' 



(r) 



Jo 



k^dk 



■ji{kr)p^xy{k). 



(79) 



Since / is an even integer, the above equation is a real number 
Once the multipoles of power spectrum p'xyik) are evaluated 
by either method described above and tabulated as a function 
ofk, we have a multipoles of the correlation function ^'xyi'') by 
a simple numerical integration of Eq. ( |79] l. Because the ver- 
tex resummation factor exponentially damps for high-A;, the 
numerical integration is stable enough. 



IV. RELATION TO PREVIOUS WORK 

A. Lagrangian resummation theory 

It is worth mentioning here the relation between the above 
formulas and previous results of Ref. ll23ll . in which the La- 
grangian resummation theory (LRT) with local Lagrangian 
bias is developed. The iPT is a superset of LRT. The results 
of Ref. 1 23] are derived from the formulas in this paper by 
restricting to the local Lagrangian bias and neglecting con- 
tributions from the primordial non-Gaussianity, although the 
way to derive the same results is apparently different. The 
definitions of Q„ and R„ functions are somehow different in 
Ref. ] |23|] from those in this paper The notational correspon- 
dences are summarized in Table fVlIIl for interested readers. 



In Ref. 112311 . the linear density field 5l and the biased den- 
sity field in Lagrangian space 5^ are related by a local relation 
^x(^) ~ ^('5l(?)) in Lagrangian configuration space. Fourier 
transforming this relation, the renormalized bias functions of 
Eq. (O in models of local Lagrangian bias reduce to scale- 
independent parameters. 



„(«) 



= (f<"*). 



(80) 



where F*"* = d"F/d6]J' is the nth derivative of the function 
F{6i,). Thus the renormalized bias functions are independent 
on wavevectors in the case of local bias, and we have {F'} = 



„(i) 



,(2) 



4" and <F") = 4^\ etc. 

It is explicitly shown that the results of Ref. [23] are exactly 
reproduced by setting X - Y, Sxx = 0, expanding the product 
(Fx ) in Rxx and adopting the replacement of variables ac- 
cording to the Table IVIIII In making such a comparison, the 
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This paper 


Ref. [23] 


Ref. [36] 


Ri(k) 


Rl(k)/PL(k) 


- 


Mk) 


R2(k)/PL(k) 


- 


R^iik) 


{F')[Ri(k) + R2(k)]/Pi^(k) 


- 


Rlik) 





- 


Qi(k) 


Qdk) 


- 


Qiik) 


Qiik) 


- 


Q3(k) 


Q^ik) - 6Q2(k) 


- 


Q4(k) 


am 


- 


Ql(k) 


{F')Q,(k) 


- 


Q^(k) 


(F')Qi{k) 


- 


Q's(k) 


{F")Qs(k) 


- 


Q^k) 


(F")Q,{k) 


- 


Q^o'(k) 


{F'?Qi{k) 


- 


eff(fc) 


(F'?Q9(k) 


- 


effw 


{F'}-Qw(k) 


- 


effw 


{F')-Qn{k) 


- 


ef/w 


{F'){F")Q,2(k) 


- 


effw 


{F"?Qn{k) 


- 


Si(k) 


- 


R2{k) 


S2(k) 


- 


2Ri(k)-R2(k) 


S^{k) 


- 


am 12 


Sf(k) 


- 


Q2(k) 



TABLE VIII: When the local Lagrangian bias is employed, and 
primordial non-Gaussianity is not considered, the expression of the 
auto power spectrum (X = Y) in this paper reproduces the result of 
Ref. 1 23]. When contributions from the primordial non-Gaussianity 
are extracted, the results of Ref. |36] are reproduced. CoiTespon- 
dences of the functions defined in this paper and those defined in 
Refs. I23l . l3ql are provided in this Table. The renormalized bias func- 
tions are constants in local bias models, and denoted by {F') = c]. 



and {F") 



^(2) 



in Ref. [23,]. 



-(l)r(l) 



product r^ Fy should be expanded up to the second-order 
terms in P]J,k) (i.e., one-loop terms). Thus, Eq. ( [3TT l is consid- 
ered as a generalization of the previous formula of Ref. ll23ll . 
Another previous formula of Ref. 02211 is a special case of 
Ref. |23[1 without biasing. As a consequence, setting c^ = 0, 
5xy = in Eq. OTT l reproduces the results of Ref. [,22,] . 



B. Scale-dependent bias and primordial non-Gaussianity 

Contributions from the primordial bispectrum, if any, are 
included in Sxy- In the cases of X - Y and X i^ Y - va, 
the relations between the primordial bispectrum and scale- 
dependent bias are akeady analyzed in Ref. [36] with gen- 
erally nonlocal Lagrangian bias. In the presence of primordial 
bispectrum, the scale-dependent bias emerges on very large 
scales ISD- The iPT generalizes the previous formula of 
the scale-dependent bias with less number of approximations. 
The previous formulas of scale-dependent bias ^ I60l - l62ll . 
which are derived in the approximation of peak-background 
split for the halo bias, are exactly reproduced as limiting cases 
of the formula derived by iPT 1 36 ] . It should be noted that the 
formula of scale-dependent bias in the framework of iPT is 



not restricted to the halo model of bias. Therefore the iPT 
provides the most general formula of the scale-dependent bias 
among previous work. The correspondence between the func- 
tions defined in Ref. \36\ and those in this paper is summa- 
rized in Table [Vlnl 

In this paper, the cross power spectrum of two differently 
biased objects, X and Y are considered in general. One can de- 
rive the scale-dependent bias of cross power spectrum Pxvik) 
as illustrated below. In the following argument, the redshift- 
space distortions are neglected for simplicity, although the 
inclusion of them is straightforward. We define the scale- 
dependent bias AbxY of cross power spectrum by 



PxYik) = [bxY(k) + AbxY(k)f Pu^ik), 



(81) 



where Pmik) is the matter power spectrum, and bxY(k) is the 
Unear bias factor of the cross power spectrum without con- 
tributions from primordial non-Gaussianity. In the lowest- 
order approximation, bxY(k) - [bx(k)bY(k)y^^, where bxik) 
and bYik) are linear bias factors of objects X and Y, respec- 
tively. When higher orders of AbxY are neglected, we have 



■^bxY = -jbxYik) 



APxYik) APUk) 



PUk) 



pm 



(82) 



where P^yi^^ '^^'^ ^m(^) ^"^ the Gaussian parts of cross 
power spectrum and the auto power spectrum of mass, respec- 
tively, and APxyik) and AP^{k) are corresponding contribu- 
tions from primordial non-Gaussianity so that the full spec- 
tra are given by PxY(k) = Pxyik) + APxY(k) and P^k) = 
Pg{k)+APUk). 

On sufficiently large scales, nonlinear gravitational evolu- 
tions are not important, and dominant contributions to the 
multi-point propagators are asymptotically given by [36|] 



^(2) 



f^^Hk) * bxik). 



r'^'ikuk2)^c'^'(kuk2). 



(83) 
(84) 



„(i) 



where bxik) = 1 + c\, ik) is the linear bias factor of object X. 



In this limit, Eq. ( |34| | reduces to 
SxYik)«bxik) [ 

Jk, 



cfikuk2)Bdk,kuk2). (85) 



In the lowest-order approximation and a large-scale limit, the 
predictions of iPT are given by 

Pgik) * Pdk), P'^yik) * bxik)bYik)Pdk), (86) 

APJk)^0, APxYik)^^[SxYik) + SYxik)], (87) 



and we have bxyik) - [bxik)bYik)] 
Substituting these equations into Eq. 



'^^ as previously noted, 
we have 



AbxYik) 



SxYik) + SYxik) 



4^jbxik)bYik)Pdk) 



(88) 



This equation gives the general formula of the scale- 
dependent bias for cross power spectra in general. 
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= ^A/Mi — ® — •vw + 



+ ^AAl^ \wv 



FIG. 4: Diagrammatic representation of the resummation scheme of 
CLPT 1 52]. The original CLPT does not include the effects of nonlo- 
cal bias, and can be easily extended to include them by applying the 
formalism of iPT and the resummation of this type of diagrams. 



FIG. 5: Ingredients of the displacement correlator. All the diagrams 
up to one-loop approximation are shown. These diagrams are re- 
summed in CLPT. 



In a case of the auto-power spectrum with X - Y, the 
above equation reduces to a known result (|3g], Abx ~ 
S xx(k)/[2bx{k)PL(k)]. Previous formulas of the scale- 
dependent bias in the approximation of peak-background split 
are reproduced in limiting cases of this result, adopting the 
renormalized bias functions c^' in the nonlocal model of 



halo bias described in Sec. IIIBI The integral of Eq. (ISSl l is 
scale-dependent according to the squeezed limit of the pri- 
mordial bispectrum, B]Jk,ki,k2) with k <k ki,k2. Thus, 
the scale-dependencies of the bias in cross power spectra are 
similar to those in auto power spectra. Amplitudes of the 
scale-dependent bias are different. When the primordial non- 
Gaussianity are actually detected, scale-dependent biases of 
cross power spectra of multiple kinds of objects would be use- 
ful to cross-check the detection. 



C. Convolution Lagrangian perturbation theory 

Recently, a further resummation method, called the con- 
volution Lagrangian perturbation theory (CLPT) [52], is pro- 
posed on the basis of LRT. The implementation of the CLPT 
actually improves the nonlinear behavior on small scales 
where the original LRT breaks down. The proposed CLPT 
is based on the LRT in which only local Lagrangian bias can 
be incorporated. 

Under the light of iPT, the resummation scheme of CLPT 
corresponds to resumming the diagrams depicted by Fig. 
The shaded ellipse with the symbol 'C represents a summa- 
tion of all the possible connected diagrams. The actual ingre- 
dients are shown in Fig.|5]up to the one-loop approximation. 
The corresponding function of this Figure is given by 

Au{k) = -L]'\k)L'^l\k)Pak) 



-I 



cPp 



•(l)^^^r(3) 



^^^^3 L\l>(k)t;>(k, p, -p)Pdp)Pdk) 
\ r Lf\kuk2)Lf(kuk2)PL(ki)PL(k2) 

^ Jk,2=k 

Jk, 



L'^^iki,k2)Bdk,ki,k2). (89) 



The indices /, j are symmetrized on RHS of the above equa- 
tion. This function is the same as Cij(k) in Ref. ||23I1 . and 



-Cij(k) in Ref. [22]. We refer to the graph of Fig. |5]and 
Eq. ( [89] ) as "displacement coiTelator" below. To the full or- 
ders, the displacement correlator Aij(k) is given by 



(¥i(kWj{k'))^ = -{Infslik + k')Aij{k). 



(90) 



The expression of Eq. d89] l is also obtained from this equa- 
tion, adopting the one-loop approximation in the perturbative 
expansion of Eq. (|9]l. 

Using the displacement correlator, the diagrams of Fig. 
can be represented by a convolution integral of the form 

Zj — r^" ■ ■ ■ ^'"^■'■' ■ ■ ■ ^>" '^'>>' ^^1^ ■ ■ ■ ^'">„(^«) 

= fd'qe-''""'' exp [-kikjAijiq)] , (91) 

where k is the wave vector of the nonlinear power spectrum 
Pxrik) to evaluate, ki...„ = ki + ■ ■ ■ k„ is the total wave vector 
that flows through the resummed part of Fig. 2] and 



Aijiq) 



J (2ny ^'^ 



(k) 



(92) 



is the displacement correlator in configuration space. The con- 
volution integral of Eq. ( 1911 1 contributes multiplicatively to the 
evaluation of the power spectrum Pxrik). 

The displacement correlator in configuration space, 
Eq. ( |92] l, is given by the full-order displacement field 'Piq) 
as 



Aijiq) ^-{nqiWjiqi))^ 



(93) 



where q = q2 - Qi- This function is denoted as Cijiq)/2 in 
Ref. 1531, and thus we have a coiTespondence, 



Cfl^^iq) = lAijiq) 



(94) 



In the CLPT, the vertex resummation factor is included in a 
function Aij{q) = Bij + Cijiq) of their notation, where Bij - 
2<T^6ij and cr^ = (|¥'p)/3. Thus we have a correspondence. 



A^I^''iq)^-cr',6u + 2Auiq). 



(95) 



The first term in the LHS corresponds to the vertex resum- 
mation in iPT, and is kept exponentiated in both original LRT 
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FIG. 6: Connected diagrams with three wavy lines up to tree-level 
approximation. These diagrams are not resummed in CLPT. 



and CLPT. The second term is kept exponentiated in CLPT 
and expanded in the original LRT formalism. 



When the Lagrangian local bias is assumed (c 



(1) 



<^'>, 



c^ - {F"},...), and the convolution resummation of Fig. |4] 
is taken into account in the iPT, the formalism of CLPT is 
exactly reproduced. When the Lagrangian nonlocal bias is al- 
lowed in the iPT with the convolution resummation, we obtain 
a natural extension of the CLPT without restricting to models 
of local Lagrangian bias. 

Extending this diagrammatic understanding of CLPT in the 
framework of iPT, it is possible to consider further convolu- 
tion resummations that are not included in the CLPT. In the 
CLPT, only connected diagrams with two wavy lines (i.e., 
Fig.|5]l are resummed. We define the three-point correlator of 
displacement, Aijk{k\,k2, kj,) where ^:i -H ^:2 H- ^3 = by the 
connected diagrams with three wavy lines as shown in Fig. |6] 
This function is given by 



D. Renormalized perturbation theory 

Recent progress in improving the standard perturbation the- 
ory (SPT) was triggered by a proposition of the renormalized 
perturbation theory (RPT) 11441 14511 . Although this theory is 
formulated in Eulerian space, there are many common fea- 
tures with iPT in which resummations in terms of Lagrangian 
picture play an important role. Below, we briefly discuss these 
common features. However, one should note that purposes 
of developing RPT and iPT are not the same. The RPT for- 
malism mainly focuses on describing nonlinear evolutions of 
density and velocity fields of matter, extrapolating the pertur- 
bation theory in Eulerian space. The iPT formalism mainly 
focuses on consistently including biasing and redshift-space 
distortions into the perturbation theory from the first principle, 
without any fitting parameter etc. The RPT (and its variants) is 
properly applicable only to unbiased matter clustering in real 
space (even though there are phenomenological approaches 
with freely fitting parameters, such as the model of Ref. |l5ll| . 
for example). Thus, the resummation methods in RPT can be 
compared only with a degraded version of iPT without biasing 
and redshift-space distortions. 



7. Propagators in high-k limit 

An important ingredient of RPT is an interpolation scheme 
between low-A: and high-A: limits of the multi-point propagator 
of mass rl^ (ki, . . ., k„) with k - k\...„. Based on the Eulerian 
picture of perturbation theory, the high-fc limit of the propaga- 
tor is analytically evaluated as [45, 46] 



{Wi{k,Wj{k2)Wj{k^))^ 



<n) 



rj{ki,...,k„) siexp 



2 ^ 



F„(ki 



,k„), (99) 



-/ (2K)^Sli(ki + k2 + ki)Aijk(ki, ^2, ^3)- (96) in the fastest growing mode of density field, where 



to the full order This three-point correlator Aijk is the same as 
—Cijk in Ref. ||23I1 and -iCijk in Ref. ll22ll . In a similar way as 
Fig.|4]and Eq. ( |9TT l, including resummations of the three-point 
correlator modifies the convolution integral of Eq. (|9Tt as 

r cfq e-''*'« exp [-kikjAijiq) + kikjkkAijkiq)] , (97) 

where k' is the total wave vector that flows through the re- 
summed part, and 

A,,i,)^l^,e-^I^A,,ik,-p,p-k). (98) 

One can similarly consider four- and higher-point convolution 
resummations, which naturally arise in two- or higher-loop 
approximations. However, it is not obvious whether or not 
progressively including such kinds of convolution resumma- 
tions actually improves the description of the strongly non- 
linear regime. Comparisons with numerical simulations are 
necessary to check. Detailed analysis of this type of exten- 
sions in the iPT is beyond the scope of this paper, and can be 
considered as an interesting subject for future work. 



o-d 



i/^^' 



PL(k'), 



(100) 



F„ is the «th-order kernel function of SPT, and k - |^i ..„|. 
Although decaying modes and the velocity sector are also in- 
cluded in the original RPT formalism [44l, M^i], we neglect 
them for our puipose of comparison between RPT and iPT. 

The multi-point propagator in the iPT has the form of 
Eqs. (|5]l and (|6]l with full orders of perturbations. In the unbi- 
ased case, X - m,we have 



rl^iki. 



<«)/ 



.,k„)^nik)r]^>iku...,k„). 



(101) 



where n{k) - (e"'* **) is the vertex resummation factor The 
Eq. ( IIOII ) should also have the same high-A; Umit as Eq. ( |99l ), 
since we are dealing with the same quantities. Although ex- 
plicitly proving this property in the framework of iPT is be- 
yond the scope of this paper, a natural expectation arises that 
the high-A: limit of the resummation factor n{k) is given by the 
exponential prefactor of Eq. (|99t , as discussed below. 

In the high-A: limit, the factor e '* ** which is averaged over 
in the resummation factor strongly oscillates as a function of 
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displacement field *?. Consequently, large values of the dis- 
placement field do not contribute to the statistical average, and 
dominant contributions come from a regime \*P\ < k^^ . In 
the high-fc limit, this condition corresponds to a weak field 
limit of the displacement field, which is well described by the 
Zel'dovich approximation, W{k) ^ {ik/k^)6h(k). Assuming 
a Gaussian initial condition, higher-order cumulants of dis- 
placement field in the Zel'dovich approximation are absent in 
Eq. ^. Since {WiWj)^ - 6ij{\V\^)ll> from rotational symme- 
try in real space, we have {{k ■ 1')^)c - k^{\V\^)/3 = k^cr^^ in 
the Zel'dovich approximation. Thus we naturally expect 



there are correspondences, 

a(k)^-ik- f^L^'^(j,)6^(j,)^-k-W'-'\ (107) 
(e-W) = (e-*^"'), (108) 

where tf'^'^ is the linear displacement field in configuration 
space at the origin. Since Eq. (llOll l holds in non-Gaussian 
initial conditions as well, the high-A: limit of iPT, Eq. ( 11021) . is 
replaced by 



n(k) = (e-*-^) * expi^-hW^ , (102) n{k) ^ (.-*«"") . exp £ ^ {{k ■ •?(»)") 



(109) 



in the high-A: limit, which agrees with the exponential prefac- 
torofEq. dllll. 

If the above expectation is correct, Eqs. d99] l and (IIOU sug- 
gests the high-A: limit of normalized propagator is given by 



which agrees with the the replacement of RPT, Eq. (1105b . 
Since only the exponential factor is replaced in Eqs. (|99] | and 
( 11021 ). the high-fc limit of Eq. ( 1103b does not change even in 
the case of non-Gaussian initial conditions. 



r'^'ik, 



, k„) X F„{ki 



,kn). 



(103) 



i.e., the high-A: limit of the normalized propagator is given by 
tree diagrams, and contributions from whole loop corrections 
are subdominant. This is a nontrivial statement, since the nor- 
malized propagator contains non-zero loop corrections in each 
order For example, taking the limit fe — » oo in Eq. (|37] l of the 
one-loop approximation, we have 



rt\k) 



58 



315 



p^dp _ 



1 + ^ -:^P^(P)^ 



27r2 



(104) 



which is apparently different from Fi - \. Actually the inte- 
gral in the RHS is logarithmically divergent for a spectrum of 
cold-dark-matter type, which has an asymptote Pi^ik) oc A:"^ 
for k ^ oo. Thus Eq. ( 1103b does not apparently hold when the 
loop corrections are truncated at any order Thus, Eq. (11031 
has a highly non-perturbative nature. This situation is natu- 
ral, because the high-fc limit of Eq. (1102b is also highly non- 
perturbative. When the equation is truncated at any order, 
a high-A: limit gives divergent terms, while the whole factor 
approaches to zero. The same is true for the high-A; limit 
in the RPT formaUsm, Eq. (|99]l. Provided that Eq. (|99]l is 
true, Eqs. ( |102b and ( |103b are the same statement because of 
Eq. ( llOlb . which is a definition of the normalized propagator. 
The above argument is readily generalized in the case of 
non-Gaussian initial conditions. In the high-A: limit of the RPT 
formalism, the exponential factor in Eq. ( |99] l is replaced by 



(105) 



exp|-lA^crd2J ^ (e-«) = exp £ £ {[a{k)r\ 



where 



2. Nonlinear interpolation I: RegPT 

In the RPT formalism, the nonlinear propagator is approxi- 
mated by analytically interpolating the behaviors in the high-A: 
limit and the low -A: limit |,45..46..49]. There are at least two 
prescriptions for the interpolation. An interpolation scheme of 
Refs. |48, 50, 51], which is called RegPT, uses an prescription 
for the multi -point propagator truncated at the A^-loop order as 



r 



RegPT 



{f„ + 6r["^ 



loop 



2-loop 



+ 6r 



in) 
W-loop 



C.T.) 



X exp 



-U^cr/\ (110) 



-('0 



where SFj^^i^^^ is the M-loop correction term of the propaga- 
tor, and C.T. is a counterterm to match the A^-loop expression 
is exact in both limits, i.e.. 



C.T. 



1,2 2 1,4 4 1 /^ 0-d^ 

-A: Vd + -k cTi + ■ ■ ■ + — 

2 8 m\ 2 



1,2 2 1 (kTd 

-k CTd + ■ ■ ■ + 

2 (A^-l)!\ 2 



5r 



Fn 



<«) 

l-loop 



+ --- + 5^W5r|;>_,,,„„p. (Ill) 



The tree-level multi-point propagators are the same as the 



kernel functions in SPT, i.e., F, 



in) 



F„. It is apparent that 



Eq. (II 10b has the correct low-A: limit. In the high-A: limit, we 
{-k^crl/2)'^FJMl according to Eq. ^. In 



have 6F 



,(n) 



/dp k ■ p 
-J-^6dpy (106) 

Comparing these equations of RPT with Eqs. ^, (|9]l of iPT, limit, F\ 



M-loop ^ ^ '^ ^ d ' 

this limit, it can be shown by induction that all the loop correc- 
tions in the first parentheses of Eq. (IllOb including the counter 
term remarkably cancel each other, leaving only the tree-level 
contribution F„. Thus Eq. ( IllOb also has the correct high-A: 



(n) 
RegPT 



F„e\p(-k^o-//2) for k ^ oo. The RegPT 
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prescription of Eq. (IllOl i can be re-expressed in a more com- 
pact form including the counterterm as 



r 



RegPT 



exp 



1 ;2 2 



^=0 



loop' 



X exp 



truncated 



-/^^ 



(112) 



-(«) 



where 5rl\^^ = F„, and [ 



trancated indicates a truncation up 
to a given order after completely expanding the exponential 
factor 

The RegPT prescription of Eq. (II 10b can be compared with 
Eq. ( llOll i in the iPT formalism. On one hand, applying a Tay- 
lor expansion of the resummation factor 77 and truncating at 
the n-loop order give the same result as the n-loop SPT. On 
the other hand, the lowest-order approximation of the resum- 
mation factor is given by Eq. (l35l l in real space, i.e.. 



77(fe)^exp(-ifeVd^|, (k- 



0), 



(113) 



which is accidentally the same as the exponential factor in the 
high-A: limit of Eq. ( |99] l. From these observations, it is now 
clear that the RegPT prescription of Eq. UIOJ is equivalent 
to evaluate the unbiased propagator /"J^ by Eq. UOIJ in the 
framework of iPT, keeping only the lowest-order term in the 
vertex resummation factor n(k) and expanding all the other 
higher-order terms from the exponent. In other words, the 
RegPT prescription is equivalent to the restricted iPT formal- 
ism where the vertex resummations are truncated at the one- 
loop level (without biasing and redshift-space distortions). 



3. Nonlinear interpolation II: MPTbreeze 

There is another scheme of interpolating the nonlinear 
propagators called MPTbreeze 114911 . which is originally em- 
ployed in the two-point propagator in the RPT formalism IHSll . 
This method is simpler than the RegPT, in a sense that calcu- 
lations of interpolated propagators require only one-loop inte- 
grals. In the MPTbreeze prescription, the interpolated propa- 
gators are given by 



-(") 



Jku- ■■ ,k„)^ F„(ku- ■ ■ , k„) exp [<'Lp(^)] . 



(114) 

where the one-loop correction term of two-point propagator 
in the growing mode is explicitly given by 



6r' 



1 -loop 



(k) 



r d'c 

J (2^ 



q Pdq) 



)3 504^3^5 



6k'' q - 19k^q^ + 5Qq^k^ 

(115) 



_ 21 V + hk^- ^2)3(2^2 ^ 7^2) i„ \kj^ 
4 \k -\- q\~ 



The notations Po(k), f(k) ff^ in Ref. l'49'l are related to 
our notations by Pi^ik) = (27ifDl{z)Po{k), ^^S'Lp(fc) = 



Dl{z)f{k) and r^^^^^^^^^ = rf/Dliz), where D^(z) is the lin- 
ear growth factor. Since 6r\ '. (k) — > -k^cr^ 12 in the high-A: 



"(1) 



1-loop^ 



limit and 5rX>^^^^{k) 



in the low-A; hmit, Eq. (11141) has 
correct limits. 

It is worth noting that the prescription of Eq. (I114t corre- 
sponds to replacing all the loop-correction terms of propaga- 
tors by 

<4oop(^l' ■ ■ ■ ' ^«) - ^ K'.oop(^)]'' Fn{k,, ..., k„). 

(116) 
Both prescriptions of MPTbreeze and RegPT give similar re- 
sults, and they agree with numerical simulations fairly well in 
the mildly nonlinear regime |49, 50]. Thus the approximation 
of Eq. ( 1116b turns out to be empirically good, although the 
physical origin of the goodness in this prescription is some- 
how unclear. 

According to Eq. ( [37] i or Eq. (|38T l, the iPT-normalized two- 
point propagator of mass is related to the function (J/"j j^^ (k) 
by 



sr?l(k) = (JfiLn (^) - :;k^(^d\ 



1-loop^ 



1 -loop ^ 



(117) 



HI) 



where dP^ j^^ is the one-loop correction term which corre- 
sponds to the integral in Eq. ( l38l l without bias, c^' - 0, or 



srtL,(k)-^Rdk) + jR2(k), 



as seen in Eq. 
have 



r 



MPTbreeze 



(118) 
HJ. Substituting Eq. ( 1117b into Eq. (1114b . we 

F„ exp [srWl^^ik)] exp I- ^k^cr/j . (119) 



Comparing this form with Eq. (1112b . the relation between the 
prescriptions of RegPT and MPTbreeze is explicit. Both pre- 
scriptions differ in the prefactor preceding to the exponential 
damping factor; a truncation scheme is employed in RegPT, 
and a simple model of the higher-loop corrections is employed 
in MPTbreeze. 



V. CONCLUSIONS 

The iPT is a unique theory of cosmological perturbations 
to predict the observable spectra of biased tracers both in real 
space and in redshift space. This theory does not have phe- 
nomenological free parameter once the bias model is fixed. In 
other words, all the uncertainties regarding biasing are packed 
into the renormalized bias functions c^ , and weakly nonlin- 
ear gravitational evolutions of spatial clustering of biased trac- 
ers are described by iPT without any ambiguity. In this way, 
the iPT separates the bias uncertainties from weakly nonlinear 
evolutions of spatial clustering. The renormalized bias func- 
tions are evaluated for a given model of bias. 

Most of physical models of bias, such as the halo bias 
and peaks bias, fall into the category of the Lagrangian 
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bias. Redshift-space distortions are simpler to describe in La- 
grangian picture than in Eulerian picture. The iPT is primarily 
based on the Lagrangian picture of perturbations, and there- 
fore effects of Lagrangian bias and redshift-space distortions 
are naturally incorporated in the framework of iPT. 

In this paper, general expressions of the one-loop power 
spectra calculated from the iPT are presented for the first time. 
The cross power spectra of differently biased objects, Pxrik), 
both in real space and in redshift space are explicitly given 
in terms of two-dimensional integrals at most up to one-loop 
order The final result in real space is given by Eq. (l47i with 
Eqs. ( I35I 1. (I42I 1. (l43l l. (I44I 1. and that in redshift space is given 
by Eq. (EB with Eqs. dST]), {5^, §0^ and (|70|. When the 
vertex resummation is not preferred, one can alternatively use 
Eq. ( |36] | instead of Eq. ( |35] |. An example of the renormalized 
bias functions is given by Eq. ( fT9] l for a simple model of halo 
bias. 

The iPT is a generalization of the method of Ref. ll23ll . 
which is applicable only to the case that the Lagrangian bias is 
local and that the initial condition is Gaussian. Although the 
derivations are quite different from each other, it is explicitly 
shown that the general iPT expression of the power spectrum 
exactly reduces to the expression of Ref. ll23ll in models of 
local Lagrangian bias and Gaussian initial condition. 

The effects of primordial non-Gaussianity are included as 
well. The consequent results are consistent with those de- 
rived by popular method of peak-background split. In fact, the 
iPT provides more accurate evaluations of the scale-dependent 
bias due to the primordial non-Gaussianity |36]. In the present 
paper, both effects of gravitational nonlinearity and primordial 
non-Gaussianity are simultaneously included in an expression 
of biased power spectrum. Thus, the most general expressions 
of power spectrum with leading-order (one-loop) nonlinearity 
and non-Gaussianity are newly obtained in this paper. 

In this paper, comparisons of the analytic expressions with 
numerical A^-body simulations are not presented. In an ac- 
companying paper ll63ll . the results in the present paper are 
used in calculating the nonlinear auto- and cross-correlation 
functions of halos and mass, and are compared with numeri- 
cal simulations, focusing on stochastic properties of bias. We 
have confirmed that effects of nonlocal bias is small in the 
weakly nonlinear regime for the Gaussian initial conditions. 
That is not surprising because nonlocality in the halo bias is 
effective on scales of the halo mass indicated by Eq. ( 1211 1; 
for example, R =^ 0.7, 1.4, 3.1, 6.6/i-'Mpc for M = 10", 
10'^, 10'"', IO'^^/z^'Mq, respectively, while one-loop perturba- 
tion theory is applicable on scales » 5-10/;"'Mpc for z < 3. 
Therefore, the predictions of iPT in Gaussian initial condi- 
tions with one-loop approximation are almost the same as 
those of LRT with Lagrangian local bias |22], which have 
been compared in detail |25] with numerical simulations of 
halos both in real space and in redshift space. The nonlocality 
of halo bias could be important beyond the one-loop approx- 
imation, which is apparently an interesting extension of the 
present work. 

In the framework of iPT, the vertex resummation is nat- 
urally defined, resulting in the resummation factor n{k) of 
Eq. (l35T l in real space or Eq. (ISTT i in redshift space. The vertex 



resummation of iPT is closely related to other resummation 
methods like RPT which are formulated in Eulerian space. 
When the vertex resummation is truncated up to one-loop or- 
der, the iPT without bias and redshift-space distortions gives 
the equivalent formalism to the RegPT, a version of RPT with 
regularized multi-point propagators. 

Beyond the vertex resummations, the scheme of CLPT 
is readily applied to the framework of iPT as discussed in 
Sec. IIVCI Further resummation scheme of convolution can be 
also considered. It might be an interesting application of iPT 
to include those type of further resummations in the presence 
of nonlocal bias and redshift-space distortions. 

Although the resummation technique has proven to be use- 
ful in the one-loop approximation, it is arguable whether the 
same is true in arbitrary orders [M]. The vertex resummation 
is not compulsory in iPT, rather it is optional. The general 
form of vertex resummation factor in iPT is given by Eq. (|6). 
When this exponential function is expanded into polynomials, 
we obtain a perturbative expression of power spectrum with- 
out resummation, which is an analogue to SPT. However, for 
evaluations of the correlation function, the exponential damp- 
ing of the resummation factor stabilizes the numerical integra- 
tions of Fourier transform, and therefore the vertex resumma- 
tion is preferred. 

The nonlocal model of halo bias 113611 explained in Sec. Ill B I 
is still primitive. There are plenty of rooms to improve the 
model of nonlocal bias in future work. The iPT provides a 
natural framework to separate tractable problems of weakly 
nonlinear evolutions of biased tracers from difficult problems 
of fully nonlinear phenomena of biasing. 



Acknowledgments 

I acknowledge support from the Ministry of Education, 
Culture, Sports, Science, and Technology, Grant-in-Aid for 
Scientific Research (C), 21540267, 2012. 



Appendix A: Diagrammatic rules 

A set of diagrammatic rules in iPT which is used in this pa- 
per is summarized in this Appendix. Full set of rules and their 
derivations are found in Ref. pi']. The relevant diagrammatic 
rules are shown in Fig. Q and |8] Physical meanings of the 
graphs are as follows: a double solid line corresponds to the 
number density field 5xik), a square box represents partial re- 
summations of dynamics and biasing, a wavy line represents 
the displacement field, a black dot represents nonlinear evo- 
lutions of the displacement field, a crossed circle represents 
the primordial spectra. The procedures for obtaining a cross 
polyspectra P^ ' ^ (ki,. . ., kj^) of different types of objects 
Xi, . . ., Xn are listed below. Auto polyspectra are obtained by 
just setting Xi - ■ ■ ■ - X^. The power spectrum is a special 
case of polyspectra with N -2. 

1 . Draw A^ square boxes with labels Xi{i - 1 , . . . A^), each 
of which has a double solid line. Label each double 
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^ n{k)c^x\ki,...,k„)k,^---k,, 



o L;"'(fci,fe2,...,fe„) 



FIG. 7: Diagrammatic rules of iPT: dynamics and biasing. 



k -k 

-* ® »— 




^ PL(fc) 



^ P<"'(fcl,fe2,...,fe„) 



FIG. 8: Diagrammatic rules of iPT: primordial spectra. 



soUd line with an outgoing wavevector that corresponds 



to an argument of the polyspectra P 






2. Consider possible ways to connect all the square boxes 
by using wavy lines, solid lines, black dots and crossed 
circles, satisfying following constraints: 

(a) An end of a wavy line should be connected to a 
square box, and the other end should be connected 
to a black dot. 

(b) An end of a solid line should be connected to a 
crossed circle, and the other end should be con- 
nected to either a square box or a black dot. 

(c) Only one wavy line can be attached to a black dot 
while arbitrary number of solid line(s) can be at- 
tached to a black dot. 



(d) A piece of graph which is connected to a single 
square box with only wavy lines or with only solid 
lines is not allowed. 

3. Label each (solid and wavy) line with a wavevector and 
its direction. The wavevectors should be conserved at 
each vertex of square box, black dot, and crossed circle. 
Label each wavy line with spatial index together with a 
wavevector 

4. Apply the diagrammatic rules of Figs.|7]and|8]to every 
distinct graphs. 

5. Integrate over wavevectors as J d^ k'. / (In)^ , where k'. 
are not determined by constraints of wavevector con- 
servation at vertices. 

6. When there are m equivalent pieces in a graph, put a 
statistical factor 1 /ml for each set of equivalent pieces. 

7. Sum up all the contributions from every distinct graphs 
up to necessary orders of perturbations. 

The rule 2. -(d) is due to partial resummations of the square 
box. For example, the left diagram of Fig.|9]is allowed. There 
is a piece of graph that is connected to a single square box 
with both wavy and solid lines. However, the right diagram 
of Fig. |9]is not allowed, because of double reasons. One is 
that the upper piece of graph is connected to a single square 
box with only wavy lines. The other is that the lower piece of 
graph with only solid lines connected to a single square box. 
Each reason itself prohibits this diagram from counted. 





allowed 



not allowed 



FIG. 9: Examples of external vertex which is allowed (left) and is 
not allowed (right) in iPT. 
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